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\Q \ We apply Bogolubov approach to QCD with two light quarks to demonstrate a 

O \ spontaneous generation of an effective interaction, leading to the Nambu - Jona- 

Lasinio model. The resulting theory contains two parameters: average low-energy 
value of Us and current light quark mass itiq. All other low-energy parameters: 
Oh. the pion decay constant, mass of the vr-meson, mass of the cr-meson and its width, 

r-| ! the constituent quark mass, the quark condensate are expressed in terms of the 

two input parameters in satisfactory correspondence to experimental data and 
chiral phenomenology. E.g. in the approximation being used we have for a^ = 
0.67 and m^ = 20MeV: U = 93MeV, m^ = 135 MeV, m„ = 492 MeV, T^ = 



Ph! 



X 



^ ■ 574 Me V, m^ = 295 MeV, <qq>= -{222 MeVf. 

PACS: 11.30.Rd, 12.38.Lg, 12.39.-x, 12.40.Yx 

1 Introduction 

The Nambu - Jona-Lasinio model PQI21E1I11 has manifested itself to be a good phenomenolog- 
ical tool for low-energy hadron physics. It is well-known, that the fundamental perturbative 
theory QCD is valid in region of large q^. In low-momenta region NJL model supplements 
the fundamental QCD. It is important, that common property of both theories consists in 
the chiral symmetry, which defines main features of low-energy hadron physics. However till 
now there was no direct derivation of NJL model from QCD. Therefore the problem to find 
a relation between parameters of NJL and those of QCD for a long time was quite actual. 
Some attempts in this direction were accompanied by inevitable introduction of additional 
parameters (see, e.g. 010]). 



In work j7] with the use of results fH^ based on the method developed by N.N. Bogol- 
ubov ini CD] an approach is proposed to obtain an effective interaction of the Nambu - 
Jona-Lasinio model, which contains no additional parameters but the QCD ones. Main 
parameters of low-energy hadron physics were calculated there [7] in chiral limit that is 
for current mass of light quarks mo = 0. In the present work we apply the approach for 
description of light zero-spin mesons with account of non-zero current quark mass. 

The structure of the work is the following. 
In Section 2 the spontaneous generation of NJL type interaction in QCD is described. 
Section 3 deals with the description of the scalar and pseudo-scalar states (mesons). 
In Section 4 spontaneous violation of the chiral invariance is demonstrated. 
In Section 5 expressions for quark condensate and for parameters of zero-spin mesons are 
presented. 
Section 6 is devoted to numerical results and discussions. 

2 Effective NJL interaction 

Now we start with QCD Lagrangian with two light quarks {u and d) with number of colours 

N = 3 

2 -. 

fe=i 

where we use the standard QCD notations. 

Let us assume that a non-local NJL interaction is spontaneously generated in this theory. 
We use Bogolubov approach jHl 1101 to check this assumption. In accordance to the approach, 
application of which to such problems are described in details in work jH] , we look for a non- 
trivial solution of a compensation equation, which is formulated on the basis of the Bogolubov 
procedure add — subtract. Namely let us rewrite the initial expression (^ in the form 

L = Lq + Lint 'i 

-i/j ip t/j tpj + — ^ ■ (^^/7^V' tpr'^'yf.ilj + V^r''757MV'V'^S57MV' j ; (2) 



- ^ • [fT Ifi'^ ^^ 7mV^ + i^T 757mV'V'^ 757m V^j • (3) 

Here ip is the isotopic doublet of quark fields, colour summation is performed inside of each 
fermion bilinear combination, -Fo^jy = dfj_Ai, — d^A^, and notation Gi ■ ipipipil) corresponds to 



non-local vertex in the momentum space 

t{27T)*GiFi{pl,p2,p3,pA)6{pl + p2 + p3 + p4); (4) 

where Fi{pl,p2,p3,p4:) is a form-factor and pi, p2, p3, p4 are incoming momenta. In 
the same way we define vertices, containing Dirac and isotopic matrices. We comment the 
composition of the vector sector, which here contain only isovector terms, in what follows. 

Let us consider expression (j2)) as the new free Lagrangian Lq, whereas expression (jH)) 
as the new interaction Lagrangian Lint- Then compensation conditions (see again jH]) will 
consist in demand of full connected four-fermion vertices, following from Lagrangian Lq, to 
be zero. This demand gives a set of non-linear equations for form-factors Fj. 

These equations according to terminology of works jHl Cni are called compensation 
equations. In a study of these equations the existence of a perturbative trivial solution 
(in our case Gi = 0) is always evident, but a non-perturbative non-trivial solution may also 
exist. Just the quest of a non-trivial solution inspires the main interest in such problems. 
It is impossible to find an exact non-trivial solution in a realistic theory, therefore the goal 
of a study is a quest of an adequate approach, the first non-perturbative approximation of 
which describes the main features of the problem. Improvement of a precision of results is 
to be achieved by corrections to the initial first approximation. 

Thus our task is to formulate the first approximation. Here the experience acquired in the 
course of performing of work ^8j is useful. Now in view of obtaining the first approximation 
we would make the following assumptions. 

1) In compensation equations we restrict ourselves by terms with loop numbers 0, 1, 2. For 
one-loop case only trivial solution exists. Two-loop terms lead to integral equations, which 
may have non-trivial solutions. So the account of two-loop terms leads to the first non-trivial 
approximation. 

2) In compensation equations we perform a procedure of linearizing over form-factor, which 
leads to linear integral equations. It means that in loop terms only one vertex contains the 
form-factor, while other vertices are considered to be point-like. In diagram form equation 
for form-factor Fi is presented in fig. 1. An accuracy of this procedure was estimated in 
work j8j to be of order of magnitude of few per cent. 

3) While evaluating diagrams with point-like vertices diverging integrals appear. Bearing in 
mind that as a result of the study we obtain form-factors decreasing at momentum infinity, 
we use an intermediate regularization by introducing UV cut-off A in the diverging integrals. 
It will be shown that results do not depend on the value of this cut-off. 

4) We use a special approximation for integrals, which is connected with transfer of a quark 
mass from its propagator to the lower limit of momentum integration. Effectively this 
leads to introduction of IR cut-off at the lower limit of integration by Euclidean momentum 
squared g^ at value m^. To justify this prescription let us consider a typical integral to be 
encountered here and perform simple evaluations. Functions which we use here depend on 



variable of the form «g^, where a is a parameter having 1/m? dimension. Now we have 



■.I2\k 



[q^ + m^f 7„2 (g'2) 

F{a g'2) dg'2 .00 p>^^ a^^a 

"^ / TTT^TI + -5 (5) 



Thus the introduction of IR cut-off at m"^ which actually consists in the following substitution 

"'^F(aq^)dq^ rF{aq^)dq^ 

-; (6) 



{q^ + ui^Y Jm^ {q^Y 

corresponds to accuracy, which is defined by parameter a m?. As we shall see, this parameter 
for our solutions does not exceed order of magnitude of few per cent. We use this tool 
throughout the present work. In doing this we keep at nominators only the leading terms in m 
expansions because taking into account of the next terms evidently means supererogation of 
accuracy. Note, that similar method was used in work |11J for description of the confinement 
in NJL model. 

5) We shall take into account only the first two terms of the 1/A^ expansion. Neglected 
terms gives contribution, which values are defined by parameter 1/(4 A^). Here additional 
factor 4 in the denominator is connected with structure of NJL interaction in Lagrangian Q. 
Indeed a trace in colour indices is always accompanied by a trace in spinor indices, which 
gives factor 4. Thus this approximation defines accuracy ~ 8 %. 

Let us formulate compensation equations taking into account all the introduced prescrip- 
tions. For free Lagrangian Lq (0) full connected four-fermion vertices are to vanish. One 
can succeed in obtaining analytic solutions for the following set of momentum variables (see 
fig. 1): left-hand legs have momenta p and — p, and right-hand legs have zero momenta. In 
particular this kinematics suits for description of zero-mass bound states. The construction 
of expressions with an arbitrary set of momenta is the problem for the subsequent approxi- 
mations. In the present work we shall use the next approximation for obtaining parameters 
of scalar and pseudo-scalar mesons. 

Now following the rules being stated above we obtain the following equation for form- 
factor Fi{p) in scalar channel 

Here integration is performed in the four-dimensional Euclidean momentum space with IR 
cut-off at mg. One- loop expressions contains terms proportional to A^ and 1 while two- loop 



terms correspond to A^^ and N. The leading terms are the same for scalar and pseudo-scalar 
cases. We perform the study with the scalar channel, because it defines spontaneous chiral 
symmetry breaking effect. Equation (jTj) evidently has trivial solution Gi = 0. Bearing in 
mind our goal to look for non-trivial solutions we divide the equation by Gi and perform 
angular integration in four- dimensional Euclidean space. From (|7j) we have 

^ ' 87r2 V ^ A2 2 2xJ 32 tt^ 

— (y^- 3/i')Fi(y) dy + - yFi{y) dy + log a; / yF^iy) dy + 

+x\ogx j Fi{y)dy+ y log y Fi{y) dy + x \^\og y + -J Fi{y) dy + 

+ ^f ^*+ (^^^ -h) [FM^y-\ [yFAvHv - 

-\ogA^(^j yF,{y)dy + xj F^{y) dy^Y, (8) 

GjNA^ / 1 GiN / 1 \ /"" r. / N ^ ^ 



^ = ml; X = p^] y = q^ ■ 



Equation (jHJ by a sequential six-fold differentiation reduces to the following differential 
equation 

rf2 / rf2 d^ f ■. f5ml ^, \\ ^Fi(x) 



dx"^ \ dx"^ V dx^ V / 4 J J X 

{Gl + 6G^G2)N 
^ 167r4 

with boundary conditions to be formulated below. 

Equation (jH)) reduces to Meijer equation [T2j. Namely with the simple substitution we 
have 



px^ 1 - VI - 64mo , 1 + VI - 64mo Pmf, 

z = — — : a = : b = : Un = . 

26 ' 4 ' 4 ' " 64 

Boundary conditions for equation (jlOj) are formulated in the same way as in works [7| , [H] . At 
first we have to choose solutions decreasing at infinity, that is combination of the following 
three solutions 



F,{z) = CiG^°(z|l, ^, ^, 0, a, 6) + G^Gf^f^z]!, ^, 6, a, ^, 0, ) + 

+ C3G^°(z|l,0,6,a, i, i); (li; 



where 



tlq 



is a Meijer function [12] with sets of upper indices a^ and of lower ones bj. In case only one 
line of parameters is written this means the presence of lower indices only, n and p in the 
case being equal to zero. 

Constants Cj are defined by boundary conditions 

O /^ lO /"OO poo POO 

-^-|/ F^{y)dy = 0■ / yF^{y)dy = 0■, / y'F,{y)dy = 0- (12) 

^TT^ 2 J^2 J.^2 J^2 

which one obtains from integral equation (jHJ by considering asymptotic behaviour of integral 
terms at infinity. These conditions and condition A = as well provide cancellation of all 
terms in equation (jH)) being proportional to A^ and log A^. Thus the result does not depend 
on a value of parameter A. By solving linear set ()12|1 . in which solution (jllll is substituted, 
we obtain the unique solution. Value of parameter mq, which is connected with initial quark 
mass, and ratio of two constants Gi we obtain from conditions Fi (yu) = 1 and 

Fi{uo) = 1; Mo 



pm^Q AT (G? + eCiGa) m^ 



26 1024 7r4 

The last line here presents the obvious condition of normalization of a form-factor on the 
mass shell. Now relations (|13|) give for A^ = 3 with the account of the first of conditions (|12|) 

uo = 1.92 ■ 10-« - 2 ■ 10-^ ; d = -| G2 . (14) 

So Gi and G2 are both defined in terms of rrio. Thus we have the unique non-trivial solution 
of the compensation equation, which contains no additional parameters. 
The form-factor now reads as (fTT|) with 

Ci = 0.28; C2 = -3.66-10^®; C3 = -7.8-10"®; (15) 

In what follows we use the notation Fi{z) for expression (fTTjl . where z is always the di- 
mensionless variable defined in eq. ()10|) . We have Fi{uq) = 1 and Fi{z) decreases with z 
increasing in the following way 

Fi{z) -^ — exp(- 3(1 - tVS) zs) + h.c. ; 

Zfi 



where D is a complex constant. It is important, that the solution exists only for positive G2 
and due to (jT^ for positive Gi as well. 

At this point we would comment the problem of accuracy of our method of taking into 
account of quark mass itiq. A possible corrections being proportional to rn^ correspond to 
dimensionless variable (jSj) where a = y/j3/8 

ami = ^/u^ ~ 10~^ ; (16) 

and so they are not significant for definition of form-factor Fi{z). 

We would also comment the composition of the vector sector. For a non-trivial solution 
with Gi 7^ we calculate one-loop terms giving contribution to equation for form-factor of 
vector terms. As a result of the first approximation we obtain just the isotopic vector terms, 
which are presented in expression (jS)). 

Let us note, that at this stage we have two possibilities: trivial solution Gj = and 
non-trivial one (jlll IT^ n3|) . We shall see below that a choice between the possibilities will 
be determined by QCD interaction. 

3 Scalar and pseudo-scalar states 

Now with the non-trivial solution of the compensation equation we arrive at an effective the- 
ory in which there are already no undesirable four-fermion terms in free Lagrangian (j2)) while 
they are evidently present in interaction Lagrangian Q. Indeed four-fermion terms in these 
two parts of the full Lagrangian differ in sign and the existence of the non-trivial solution of 
compensation equation for Lagrangian Q means non-existence of the would be analogous 
equation, formulated for signs of four-fermion terms in interaction Lagrangian Q.^ 

So provided the non-trivial solution is realized the compensated terms go out from La- 
grangian (0) and we obtain the following Lagrangian 

Q / _ _ _ \ 

- -Y ■ (V'^SmV' V'^SmV' + V'^SsTmV'^^SsTm^J ; (17) 

where Gi, G2 are defined by relations (fT^ and form-factor Fi is defined by eqs. (fTTllT^j) . 

Here we have to comment the meaning of the strong coupling constant g,^. It is well- 
known that g'^/4:7T = as{q'^) is the running constant depending on the momentum variable. 
We need this constant in the low-momenta region. However the perturbation theory in QCD 
does not work for small q^. We assume that in this region as(g^) may be approximated by 

^In other words the fact, that sum of a series ^ G"^an = for some value G, by no means leads to a 
conclusion, that sum of the same series with G —* ~ G vanishes as well. 



its average value ag. This assumption is very close to conception of frozen strong coupling 
at low momenta JSj • The consideration of average low-momenta value of a^ P^ El 1^ 1^ 
leads to definition of a possible range of values of a^ from 0.45 up to 0.75. So in what follows 
we use constant a^ which is assumed to fit this interval of possible values. 

Thus, bound state problems in the present approach are formulated starting from La- 
grangian (fTTj) . 

Let us write down Bethe-Salpeter equation for a state in zero-spin (scalar and pseudo- 
scalar) channel in the same approximation as was used in equation (|7j). Let us begin with 
massless states. The definitions of momenta are the same as in eq.(|7|) 



GiN r m{q^)dq {Gl + QGi G2)N 
2 7r4 J ^2 + ^^ 



^(p') = ^ / ^^^^ + ' '\7_r' X 



oa2 , / n2, i.P-(lf 3 2 ^^ \^{q^)dq 

2A +(p-g) log -iP-(l) -777 T? 5 • (^^ 



A2 2'"^ ^' 2{p-qfJ q 

Here m is a quark mass, which in general may not coincide with ttiq. We define the value of 
m after considering the spontaneous breaking of the chiral symmetry. 

After angular integrations we obtain the one- dimensional equation similar to eq.(|H|) 

^ , , Gi AT p {Gl + QGiG2)N [I r . ^ ^ 4.^, ., 

27r^ y„2 32 vr* \bx J^2 

o rx px px poo 

+ 7; y'^{y)dy+ log X y^{y)dy + x\ogx I '^{y)dy+ y\ogy'^{y)dy + 

^ J m? J rrfi J m? Jx 

+ xj (log y +-)^{y)dy + / -^!il dy + (2A^ - -x) J ^ ^{y) dy - 



/ y^!{y)dy + x '^{y)dy 

J m? J m? 



The corresponding differential equation for \E'(a;) is almost the same, as the previous one 
with one essential difference. Namely the sign afore j3 is opposite. 

dz )\ dz ) dz\ dz 2)\ dz 2)\dz ) 2^ ^' 



2 -1 + VI + 64n - -1- V1 + 64U I3m^ 

z = x- -a= ^ ; h= ^ ; « = -^ • (19) 

In this case we have the following solution decreasing at infinity 

1.1 .\ „. „,n / , . 1 1 



<^{z) =CrG^°(z|l, -,0, -,a,6) +C*G^°(^|1, -, -,0,a,6) + 

+ ClG'^[z\l, a, 6, 1 ^, 0) + ClGf,[z\\, a, 6, 1, ^, o) . (20) 



z 



Constants C* are defined by the following conditions 

/ ^ ^ = 0; / '^{z)dz = 0; / yG <if (z) dz = ; ^(m) = 1. (21) 

where u = j3rn^ /2^. Let us remind, that boundary conditions (PT|) guarantee cancellation of 
terms in equation (|T8|) containing cut-off A. Performing integrations in expressions (J2T| . we 
have the following set of equations 

ClG^[u\l, i, 0, i, a, h) + C*Gg°(«|l, 1, ^, 0, a, fe) + 

+ C*Gf,[u\l, a, 6, ^, ^, O) + C:G30(^^|1^ ^^ ^^ 1, 1, o) = 1; 

GiGoel^nl^,!, ^,0, 2+a, 2+^j + ^^ U(l/2)r(l/2 - a)r(l/2 - 6) 

— G^-j^U 13/2,1,1,1/2,0, l/2+a,l/2+fej j ~ ^3 ^06 (^ "" I g ' 2 ~^ ^' 2 ~^ 2 ' "*"' ^/ ~ 

-QGg°(«|l, i + a, i + 6, |i,0) =0; 

^ Vr(-i/2)r(-a)r(-6) ^^r i2,3/2,i,3/2,o,i+a,i+6; j 

-G;G'^(^u\2, 11 0, 1 + a, 1 + 5) - C3*Gg°(«|2, 1 + a, 1 + 6, | | o) - 
-GlGf,[u\^-A + dA + \'2,\.^) = 0; (22) 

r(5/2) 



'^rG'?7(w 11,2,1, 2,0, |+a,|+fe) + '^2*(^ 



r(-l/2)r(-l/2 - a)r(-l/2 - b) 



■-^171 " 15/2,2, 2, 3/2,0, 3/2+a,3/2+fe/ j ^3 "-^irl " l5/2, 3/2+a, 3/2+6, 2, 2, 3/2, 



— c* n^^l 11 1^ - I — n 

*-"4 "-^17^^ " 1 2, 3/2+a, 3/2+6,5/2, 2,3/2, Oy ^ ' 



For a given value of u these conditions (j^ uniquely define four coefficients G*. The result, 
that equation (fTHj) has unique solution, which satisfies all boundary conditions, corresponds 
to existence of a zero-mass state in the same approximation as is used for compensation 
equation ((7j). This is quite natural due to Bogolubov-Goldstone theorem p inH ITT]. 

However we have to take into account chromodynamic interaction as well as an interaction 
of these mesons (0 and tTq) with quarks. Indeed we have just shown the existence of this 
states and so the following effective meson-quark interaction is to exist 



(pipip + tiraip'y^Tai'y, (23) 



where g is defined by normalization condition of zero-spin states 



Tlie form-factor of interaction (j^Hj) for our standard quark momenta prescription {p, — p) is 
Bethe-Salpeter wave function defined by eqs. ((201 1221) • The account of contributions of meson- 
quark interaction was considered in the framework of the Nambu - Jona-Lasinio model e.g. 
in works fHlE] and was shown to be corresponding to the next order of the 1/A^ expansion. 
Let us calculate a mass correction term due to these contributions. For the purpose let 
us take into account terms of the first order in P^, where P is the momentum of a scalar 
(and pseudo-scalar) meson and one-loop terms being due to quark-gluon QCD interaction 
and quark-meson vertices. Note that for the last loops we use massless meson exchange. 
We define momenta of left-hand legs in fig. 2 to be p + P/2 and —p + P/2 and obtain the 
following equation 



*.(/)^%?/^(i 



3P2 (gp)2. (Gj + 6G^G2)N 

H — '/ — ^TTTT I + ^tt: — c ^ 



4g2 (q2)2 J 32 7r6 

2A' + ip~griog^^^-l-ip-gY-—^—^]{l-:^+^ii^]x 



A2 2'"^ ^' 2{p-q)V\ 4g2 (g2)2 

Jilpiq^)dq ^ gl g^ \ f ^ p{q^) dq 

g2 ^ Utt^^SttV J q2(q-p)2- ^ > 

In the course of QCD term calculation we use transverse Landau gauge ^ . Let us multiply 
equation (j^ by \l/p(p^)/p^ at P = and integrate by p. Due to equation (|THj) be satisfied 
we have 

P2 fm[q^fdq ^ f^^ JJ^\ f ^{p')dp r ^{q')dq ^ ^ ^^q) 



2 J (g2)2 V47r4 SvrV J p'^ J q^{q-pf 

After angular integration we get 

2 8 y„2 \x J^2 J^ y 

'29s +9^ r ^{z)dz [' m{t)dt {2gl + g^)h P x^ (3y^ 



2^ Ju z J^ ^ 2v^ ' 64 ' 64 

Integral inside of I^ with account of boundary conditions ()21|) reads 



(27) 



r-^--c'M^^wi^^,..\.-^ 



3 11 1 7\ ^* ^^n / , 3 1 It 1 

|0,l,-,-,- + a,- + 6)+C3*G-(z|-,- + .,- + M,2^ 
+ Q*G|^°f;2|l, ^ + a, ^ + 6, 0, ^, ^); (2J 



2 '2 2 2 



^In the approximation used the transverse gauge leads to absence of renormalization of both vertex and 
spinor field. 

10 



and after substitution of relation (j2oj) into integral J5 it is to be calculated numerically. 
Note that while evaluating integral ()28|) we use the following relation, which is presented in 
work [3 

G?K-lt.«...') = r(i-/r'(i-!!r(i-i) " <'^'>' =• '• «■ "■ ') ' f^^' 

Integral J5 turns to be positive, so the mass squared of scalar and pseudo-scalar mesons 
is shifted to negative value 

2 f(^s , ^^ \ 8/5 

4 Spontaneous breaking of the chiral symmetry 

The negative value of m^ (J3U|) means instability of the vacuum. Therefore we have to consider 
an effective potential depending on scalar field (p. In doing this we need an expression for 
mass operator of the quark E(p^). The Schwinger-Dyson equation defining this function in 
our approximation reads as follows 

GiN f ^{q^)dq {Gl + QGi G2)N 
2n^ J ^2 + ^^ 



W) =mo+^^ I ^^^H^ + ^"^Vo_r^^^^ X 



+ (^ + JL\ f ^(g')^g (31) 

Un^ Stc^J J q^iq-pf ^ ' 

The first approximation corresponds to mo = Qs = 9 = ^- Then equation (|^ exactly 
coincides with equation ()18p for Bethe-Salpeter wave function \l/(p2) ()20|) . Similar situation 
takes place in standard NJL model 0. For non-zero mo we have without gluon and meson 
corrections 

S(x) = mo + (m, — ?7io) ^E'(x) ; S(— m^) = m. (32) 

Emphasize that approximate solution (j32p of equation (j31|) exists for any value of m. For 
definition of m one has to turn to the spontaneous breaking of the chiral symmetry. 

Let us write down the effective potential which defines a possibility of the symmetry 
breaking. We look for terms proportional to 0" for n = 1, 2, 3, 4. The term with n = 2 is 
evidently defined by (jHUI) . For terms with n = 3, 4 we take quark-loop diagrams with three 
and four scalar legs respectfully and as a result we have the following effective potential 

4/^ /. /%^'/^ / 1 , a, ^ Js^ / 1 , a, ^ he 






11 



Here 

The connection between terms with n = 1 and n = 2 is obtained from the fact, that the 
tadpole term due to eq.(|H^ gives just the same contribution as the two-loop one up to factor 
{m — mo)/g- The contribution to the tadpole term being proportional to mo is zero due to 
boundary conditions ((211) • 

As for one- loop terms with n > 5, they all converge with point-like vertices. In this case 
they can be calculated and summed up to give the following additional term 

^y = ^^TTT^a-o'iogii-ei + t/^ - ^^ + 14 - ^)-^ (35) 

VlbTT"' IbTT^ 62 TT^ 48 TT^ 192 vr^/ 

which evidently does not destroy stability conditions and turns to influence results quite 
insignificantly. Thus we neglect it. 

We look for a minimum of potential (jH^ that is for a solution of the following equation 

Constituent quark mass is expressed through the vacuum expectation value of scalar field 

m = mo + grj; t] = < (p > . (37) 

Bearing in mind definitions (fT^OUl) of parameters mq and u, we come to the conclusion, that 
the position of minimum ^o has to be the following 

Thus from relations (j^ IHHl ESI EH) we obtain the following expression for a^ 

Here all integrals are functions of u and so relation ()39p defines function as{u). 

Now it is the proper place to comment the problem of stability. From the very beginning 
we have two solutions: the trivial one Gi = G2 = 0, m, = mo and the non-trivial one, which 
in details is presented above. We get convinced that the non-trivial solution corresponds to 
the minimal negative value of effective potential (j33p while the trivial solution corresponds 
to its value zero. So we are to conclude, that just the non-trivial solution is stable and thus 
the non-trivial solution is to describe the observable physical quantities. 



12 



We apply quark mass operator ()32|) to obtain also the expression for pion decay constant 
/tt. Considering one-loop quark diagram for decay amplitude of process vr"*" — >• /i"*" u^ we have 



A = f^ / ((m-mo)^(i/)' + mo^(y))— = 

= j^({jn-mQ)h+mohy, h = ^^ • (40) 

Provided either ttiq = ot I2 = h we get with account of normalization condition (j^^ just 
the original Goldberger - Treiman relation m = g f.^. We use full relation pn)) . However let 
us note, that values of the two integrals are close I2 — h and the simple original relation 
works with sufficient accuracy^. 

5 Pion mass and quark condensate 

In relations (jHH| En|l we have used approximation ()H2|1 for the quark mass operator. For 
calculation of the pion mass and the quark condensate we need the next approximation for 
mass operator. In view of this we reformulate equation (|HT|) for the following function 

$(/) = ^^^^1^^ ; (41) 

m — niQ 

and the first approximation for $ is just \1/. Then we introduce (j^T| into (J3T|l to obtain 

, 2^ ^ G^ r {{m-mo)<l>{q^) + m,)dq (G^ + 6GiG'2)iV 
^^ ^ 27r4 J {m-mo)q^ 32 7r6 ^ 



A^ 2 2{p — q)'^/ {m — nio) q 

, (jI_ , _9^\ f (mp + (m-mo)$(g^))rfg 
Uvr^^STrVy (m-mo)g2(g-p)2 " ^ ^ 

Now we subtract equation ()25|1 from equation ()42|1 and obtain the following relation 

2^ ^ G,N ( n{m-mo)D{q') + rap) dq ^ C ^H jq') dq f 3 P' (gP)^^ ^ 
2 7r4 Vi (m-mo)g2 7 g^ V4g2 {q^Y ) ) 

((m-mo)£'(g^)+mo)rfg f '${q'^)dq /3P^ _ (gP)S 

[m — mo) q^ J q^ \Aq^ [q^Y ) 

^Additional contributions to f-n being due to the P-exponent, which maintain the electroweak gauge 
invariance, are shown to cancel in the chiral limit 20. . 
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,(^,^\ f {{m~mo)D{q'^) + mo)dq 
\An^ 8 7T^) J (m-mo)g2(g-p)2 " ^ ' 

The analogous procedure is applied in standard NJL model 1^, [3| while proving the Gell- 
Mann, Oaks and Renner theorem f^. Then we again multiply ()43|) by \l/(p^)/p^ and in- 
tegrate over dp. Due to equation ()25|) be satisfied only terms being proportional either to 
P^ = — m^ or to mo do not cancel and finally we have the following relation for mass of the 
TT-meson 

2 vr (m — mo) -(2 V w V Stt/ 



'log 



[^^'^{z)dz = ClG^{u\\,\,Q,\ + a,\ + h,Q) 



-C;Gf,{u\^-,Q,Q^-^- + a^- + h) -QG^O(«||i + a, 1 + 5,1,0,0) + 

+ C:G^°(«|0, i + a,l + 6, |1,0). 

We see that pion mass squared is proportional to rriQ in accordance to the result of well- 
known work j2Tj. Note, that contributions being proportional to ttiq arising from the first 
two terms of equation (jlH|) are summed to overall zero due to consequences of boundary 
conditions (J22I)- 

From equation (PT|) we obtain also the next approximation for $, which leads to a non- 
zero value of the quark condensate 

AN f ^iq)-mo ^ 
[2tt)^ J g^ + m^ 
_ iV(m-mo) /a.^^\ r dz ^ r ^<{t)dt ^ r^ ^>{t)dU _ 



71^(3 \n 871"^/ J^ ^/z \J^ v^ Jz t 

N{m-mo) /as 9^ \ [ T *W logtrft ^ ^ [^-<^{t)dt 



After evaluating the integrals we have 

/ q"^ \ 3 m^ (m — rrio) / ^^ ( ^-xn, ,3 1 1 1 ,, 

<..>--(«. + |^)^^^(-Cr(G30(.|--,0,0,- + a,- + 6) 



-2Gg°(«|l, i,i,0, i + a, ^ + 6)) + 



C*(G^°(n||o,0,^,^ + a, i + 6) + 2G^°(«|l, l,i, i, i + a, i + 6)) + (46) 
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-1-1 OT 11 "IT 

-CI {Gfei^^\- + a,- + b,0, -, -, 0) + 2Gll{u\l,- + a,- + b, 1, -, - 

Scalar field corresponds to the cr-meson. To calculate mass of the a-nieson we use 
relation for difference of a and vr masses squared following from one-loop diagram 

2 g'N fm^nql, ,.^, 

^'^-^^ = ^^1 ig2 + ^2y dq. (47) 

Following again our rules and using expression ()32|1 we have for mass of the a-meson 

ml = ml ^ — (ml h + 2mo(m - mo) h + {m - m^f J4J ; (48) 

The 0" vr TT vertex gives according to triangle one-loop diagram the following coupling constant 

(niQ J3 + (m - mo) I^ ; (49) 






7r2 



and the cr-meson width reads 



'--- 'ml -Ami. (50) 



16 71 ml 

6 Numerical results and discussion 

Now we have expressions for all quantities under study. Then we proceed as follows. 

1) We calculate function a^ (IHUj) depending on parameter u ((201) and get convinced, that the 
interesting range of a^ corresponds to u varying in the following region 

0.0005 < u < 0.0015. (51) 

In doing this we use parameter uq = 2 10^^ according to relation (jl4|) and calculate constants 
C*, z = 1, 2, 3, 4 from boundary conditions ()22|) thus defining "^{z). Having "^{z) we calculate 
integrals Ij, j = 1, 2, 3, 4, 5. 

2) We fix value /^ = 93 MeV. 

3) Then for given u in range ()5H) from ()40|) we obtain constituent quark mass m. 

4) Having m and a^ we calculate m^ from (j44j) . 

For u in range (j51|) m.^^ varies insignificantly between 134 MeV and 135 MeV with maxi- 
mal value 134.8 MeV" at m = 0.0009, that corresponds to a, = 0.673 and mo = 20.27 MeV. 
Considering this value of ?7i,r to be the most suitable, we present a set of calculated parame- 
ters for this conditions including quark condensate ()46|) and parameters of the cr-meson 
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p]|) as well 

a, = 0.673; nio = 20.3 Mel/; 

m^ = 135 MeV; f^ = 93MeV; m^ = 492 MeV; T„ = 574 Mel/ 

m = 295 MeV ; <qq>= - (222 MeVf ; (52) 

G^ = --r ; q = 3.16. 

^ (244Me\/)2' ^ 

The upper line here is our input, while all other quantities are calculated from these two 
fundamental parameters. We present here also values of four-fermion constant Gi and of 
meson-quark coupling g. Note, that for these calculations uncertainties due to our method 
of infrared cut-off are defined not by eq. (fT^ but by the following quantity 



u 



VO.0009 = 310" 



that is the accuracy of numbers ()52|) is evidently not better than 3%. There are also other 
sources of uncertainties and so we may estimate the overall accuracy to be of order of 10%. 
The main contribution to this estimate is provided by the next orders of 1/A^ expansion, 
according to the discussion in the Section 2. 

Bearing in mind the last remarks, we may consider the correspondence of our results for 
m^, /tt and < g g > to existing data being quite satisfactory. Value for constituent quark mass 
is also consistent. As for parameters of the a-meson, experimental data according to [22j give 
a wide range for their possible values. However let us note, that recent determinations [23|l2^ 
of the (T-meson parameters give more definite results. They are respectfully the following 

m^ = 470 ± 30 , r^ = 590 ± 40 ; m^ = 541 ± 39 , T^ = 504 ± 80 ; (53) 

that agrees with the present calculations. There is a recent analysis of the vr — vr data with 
light a [22] 5 which also agrees with a parameters (j^ . 

To conclude we would like to emphasize that the present approach for the first time 
permits to determine parameters of effective interaction inherent to the Nambu - Jona- 
Lasinio model in terms of parameters of the fundamental QCD. The optimal value of a^ = 
0.67 in (|52p is quite reasonable from the point of view of the existing knowledge on its low- 
momenta behaviour (see again f^). As for value of current quark mass mo ~ 20MeV^, 
it seems to be rather larger than usual values mo{2GeV) ~ 4 — 8 MeV. To comment 
the situation let us note, that firstly the low value of mo being mentioned corresponds to 
perturbative region and the problem how this running parameter varies while yg^ moves 
to low energy region deserves a special study. In considering of the running rriQ we have to 
take into account the effective NJL interaction as well. Secondly, the lattice studies give as 
a rule rather high values for mo, e.g. in work ^26.; values of mo corresponds just to few tens 
of MeV. The smaller values are to be obtained in the continuous limit, which till now is 
performed only by an extrapolation procedure. 
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So we may state that the aim of the work is achieved. We have begun with the demon- 
stration of the non-trivial solution of the compensation equation. The appearance of scalar 
and pseudo-scalar excitations (mesons) in the same approximation is a consequence of its 
existence. The account of QCD interaction and of meson-quark interaction leads to the shift 
of their masses squared to the negative region, i.e. to the appearance of tachyons, which are 
necessary for scalar condensate to arise. As a result we obtain the standard scheme leading 
to the spontaneous chiral symmetry breaking. Subsequent approximations of the approach 
are related to values of the quark condensate and of the pion mass. 

We have shown that the application of the method of works [7| |H], which is based on 
Bogolubov approach, to the low-energy region of hadron physics leads to quite reasonable 
results. Let us once more emphasize that we have no additional parameters but those 
entering in the low-energy QCD: as and rriQ. Thus we derive effective interaction of NJL 
type from the fundamental QCD. Basing on the results of the work we would make two 
essential conclusions. 

Firstly, a subsequent development of the present approach to the hadron physics quite 
deserves attention. In particular it is advisable to apply the approach to calculation of 
parameters of vector mesons p, u;, Ai, to consider hadrons containing s-quark, to take into 
account the vr — Ai-mixing, to study diquarks etc.. These problems comprise subjects for 
forthcoming studies. 

Secondly, the positive result of applicability test with Nambu - Jona-Lasinio model, 
being taken as an example, allows to hope for successful application of the approach to other 
problems. In particular we mean the problem of a dynamical breaking of the electroweak 
symmetry. A qualitative discussion of possible variants in this direction is presented e.g. in 
work 127]. 
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Figure captions 



Fig. 1. Diagram representation of tlie compensation equation. Black spot corresponds to 
four-fermion vertex with a form-factor. Simple point corresponds to a point-like vertex. 

Fig. 2. Diagram representation of Bethe-Salpeter equation for scalar and pseudo-scalar 
bound state. A meson corresponds to a double line. Dotted line corresponds to the gluon. 
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